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SUMMARY 


Classical drag equations in a modern version have been used 
to calculate the influence of tail modifications on the drag of a 
standard class sailplane. The profile drag which depends on the 
Reynolds number is included in the calculations. Minimum drag is 
compared with real drag for two lift coefficients. 

Some results have no clear tendency but low tail area and 
relatively low tail aspect ratio give some advantages. Optimum 
and real lift ratios between wing and tail plane are compared for 
the original sailplane. 


INTRODUCTION 


Since the energy crisis in 1973 there is a lot of interest in 
reducing the trimmed drag of airplanes (Refs. 1-5). One contribu- 
tion to the trimmed drag is the wing/tail interference drag. This 
interference drag had been interpreted as a component of the t^il 
lift vector due to local downwash angle at the tail position. Sachs 
(Refs. 6, 7, 8, 9) has shown that this interpretation is not cor- 
rect. The exact method is to calculate the interference drag with 
the aid of the downwash angle at downstream infinity. This new 
explanation corresponds to the well-known biplane theory of Prandtl 
and Munk (Ref. 10) which was also used in some new papers (Refs. 4, 
11 , 12). This theory in the modern version was used in this paper 
to show the relation between optimum and real load distributions 
between wing and tail. 

Another purpose of this paper is to show the influence of 
tail plane design on total drag. It is important that the 
addition of the Reynolds number dependent profile drag has great 
influence on the optimum design. 

All calculations are performed for a standard class 
sailplane . 
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SYMBOLS 


Values are given in SI units. 


AR 

aspect ratio 

b,b t 

wing span, tail span, m 

C D 

drag coefficient 

* 

C D 

drag coefficient (tail profile 
drag) 


drag + total induced 


C 


L 


lift coefficient 


C m pitching moment coefficient 

c mean aerodynamic chord, m 


D 

g 


drag, N 

. 2 
acceleration of gravity, m/s 


h ' h wb' h t 


distance in chord lengths from leading edge of wing 
to c.g., wing-body aerodynamic center, tail 


b wb' b t induced drag factor for wing-body, tail 

L lift, N 


m 

q 

Re 

s,s t 

e 

£ 


u 

a 


mass, kg 

2 

dynamic pressure, N/m 
Reynolds number 
wing and tail area, 
downwash angle 

downwash angle at downstream infinity 

downwash factor, e = e* k , C T 

wb L wb 

span ratio b^/b 
interference factor 
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Subscripts: 

min 

o 

opt 

t 

wb 


minimum 
zero lift 
optimum 
tail 

wing-body 


BASIC RELATIONS 


The following fundamental relations were used, assuming that 
the_aerodynamics are linear and that the dynamic pressure ratio 
is q t /q = 1: 


c + c 

L . ^L. S 

wb t 


: = C + C T (h - h , ) - C T (h, - h) = 0 

m mo . L , wb L, S t 

wb wb t 


C = C + C T (h - h , ) - C T (h, - h , ) = 0 

m mo , L wb L. S t wb 

wb t 


2 t 2 

’ = + k , CZ + -=r (r + k, C T + e C T ) 

D Do , wb L , S Do. t L, » L. 

wb wb t t t 


( 1 ) 

(2) 

(2a) 

(3) 


The last term within the round brackets of the drag equation (3) 
is the wing/tail interference drag. The derivation was given by 
Sachs and shall not be repeated here. This interference drag is 
the product of the tail lift and the downwash angle at downstream 
infinity and corresponds to the one given by Prandtl but the new 
expression is much easier to use in calculations. The downwash 
angle at downstream infinity may be expressed as (Ref. 13) 


' E * k „b C L ^ <‘* = 0 - 3) < 4 > 

wb 

The downwash factor e* = 1 corresponds to a rectangular, e* = 2 
to an elliptic and e* = 3 to a parabolic spanwise lift 
distribution of the wing. 
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Influence of Reynolds number 

Wortmann (Ref. 14) has designed and measured a lot of 
excellent profiles for sailplanes. The wind tunnel test results 
are published in the "Stuttgarter Prof ilkatalog" . Figure 1 shows 
some test results demonstrating the influence of Reynolds number 
on profile drag for several profiles. The solid lines are 
according to the relations used in this paper for calculating 
the influence of Reynolds number: 


'Do 


wb 


'Do, 


0.009 

Re 0 ; 3 

wb 

0.007 
0.3 


Re 


Re in millions 


(5) 


These relations are only valid for the above mentioned profiles 
and for a special Reynolds number range. 


Real lift ratio 


( 2 ) 


For balance in equilibrium flight (C m = 0) equations (1) and 
can be solved for the lift ratio between tail and wing 




= const 
= 0 



'mo 


wb 


+ (h - 


h wb> 


( 6 ) 


mo 


wb 


+ (h " h wb> 


- (h, - 


h wb } 


This lift ratio depends on fixed quantities and on the parameters , 
total lift coefficient and c.g. position. It is possible to elimi- 
nate the c.g. position by considering the stability requirement. 
The static margin may be expressed as 


9C 

3C 


m 

L 


(h - 


h wb> 


(C L V S t 

-cr^ -r (h t - 


h . ) (1 - ~) 
wb 3 a 


Solving equation (7) for (h - h and combining this with 
equation (6) gives the following equation: 


(7) 
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wb 


3 C 
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wb 


9C t 


(C L > S 

-CT-^IT (h_1 W (1 "lf ) 


J a 


C T =const 

C ~=0 
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mo 


wb 


3C 
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(C L > S. 

at t 


( 8 ) 


3C, 


< h - h wb»' 1 -|f>-<V h wb> 


Regarding equation (8) it is possible to say that the real load 
distribution between tail and wing depends on several fixed values 
and on the parameters, stability margin and total lift coefficient. 


Optimum load Distribution 

Prandtl has published the optimum tail/wing lift ratio in 
his biplane theory (Ref. 10) : 


wb opt 



with vi 



(9) 


The interference factor a has to be taken out of diagrams (Ref. 6). 
Therefore it is easier to use the following relations. Comparing 
Prandtl 1 s equation for the interference drag 


D 


2 a L t L wb 

Int " - b. b 

it q t 


( 10 ) 


with equation (3) it is possible to rewrite equation (9) as 

- * 


L t 

{ fr~ ) 

wb opt 


1 - 


(— ) - — 

V 2 


( 11 ) 


The optimum tail lift may be either positive or negative 
depending on the downwash factor e*. Only an elliptic spanwise 
lift distribution over the wing requires zero tail lift. The 
combination of equations (6) and (11) gives the optimum c.g. 
position 


(h - h 


wb^ opt 


mo 
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wb 
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and the optimum stability margin 

SC. 


m. 


( — — ) 


9C 


L opt 


" (h ' h wb»opt 


<C L * S 
at t 


(h, - 


h wb> (1 


l£) 

9a ' 


(13) 


a 


The minimum induced drag was also given by Prandtl and may be 
rewritten in the following equation 


’ = v n 

'Di . wb 

min 


1 - 


(1 - 


-) 


1 + 


<^> 2 - 

b t 


(14) 


The term in the brackets demonstrates the decreasing induced drag 
of the wing/tail combination compared with the wing alone. 


NUMERICAL RESULTS 


Optimum c.g. position 

All calculations are performed with the data set of a typical 
standard class sailplane; see table 1. 

The pilot of this airplane wants to know whether he can 
reach good performance by choosing the correct c.g. position.. 

In figure 2 the equations (8) and (11) are evaluated. It is easy 
to see that the optimum load distribution depends on the 
downwash factor e* but not on the lift coefficient C^ while the 
real lift ratio has inversed dependencies. The downwash factor 
is normally not known exactly; therefore it is not easy to reach 
exact conclusions from figure 2. A more accurate way is to eval- 
uate the equation (13); see figure 3. It is easy to see that 
the optimum stability margin for lift coefficients C L > 0.4 is 
obtained by a normal (stable) c.g. position. This statement is 
valid for downwash factors between 1 and 2. Good performance for 
a wide range of lift coefficients will therefore be obtained by 
choosing a medium or forward c.g. position. Only high speed flight 
requires aft c.g. positions. 


Tail modifications 

While the pilot of a sailplane is interested in the optimum 
c.g. position, the sailplane designer is interested in the tail 
design to meet stability requirements and to achieve good per- 
formance. For a fixed wing geometry it is possible to vary two 
main parameters: tail area and tail span. It is important not to 
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evaluate the parameter variations only by regarding the total 
induced drag because the tail profile drag will also change. 
Therefore the criterion is the sum of both the total induced drag 
and the tail profile drag (equation 5) 



The induced drag is calculated as in equation (3) with the lift 
coefficients obtained from the real load distribution (equation 8) . 
The stability margin is assumed to be = -0.15, and the 

downwash factor is assumed to be e* = 1. (The calculations were 
also performed for e* = 2; the tendencies correspond to e* = 1.) 

The real drag coefficient is compared with the theoretical 
minimum drag coefficient using equation (14) for the minimum 
induced drag. 

In figure 4 this drag coefficient is plotted versus the tail 
span with constant tail m.a.c. for two lift coefficients. Param- 
eter in this diagram is the tail profile drag at Re = 10^. The 
minimum induced drag (Cj) 0 . = 0) is decreasing with increasing tail 

span but the higher profile drag due to increasing tail area is 
predominant. It is suitable to design the tail with the minimum 
possible area to satisfy stability requirements. 

Assuming a minimum tail area of 1 m^, another question is, 
what span or what AR^ is optimum? There are two effects: 

- With increasing tail span the minimum induced drag of the 
complete sailplane is decreasing 

- Increasing tail span means decreasing Reynolds number resulting 
in increasing tail profile drag. 

In figure 5 it is shown that the superposition of these two 
effects results in no clear tendency. With increasing profile drag 
the Reynolds number effect becomes predominant. Assuming a tail 
profile drag coefficient of Cp 0 t = 0.01 at Re^- = 10^ a reduction 
of tail span from 2.4 m (original value) to perhaps 2.0 m will 
give some little advantages. 

The differences between minimum and real drag coefficients 
are generally small; only low total lift coefficients (C L = 0.2) 
require high tail downloads (see figure 2) resulting in greater 
differences . 


Wing and tail modifications 

The wing geometry is included in the variations. The only 
restrictions are now a wing span of 15 m and a total area (wing + 
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tail) of 11 m 2 . To simplify the calculations it is assumed that 
wing and tail have the same profile drag at Re = 10 6 . It is 
suitable to regard the drag or the drag areas directly rather 
than the drag coefficient: 


E° . c s 

5 C DO 
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1 + (£-) - e * 
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(17) 


Figure 6 shows the total profile drag area (equation 16) plotted 
versus the ratio of tail area to total area. For normal tail 
aspect ratios between 3 and 6 the total profile drag will 
increase with increasing tail area due to decreasing medium 
Reynolds number of the total area. 

The combination of equations (16) and (17) gives the minimum 
total drag which is plotted versus the ratio of tail area 

to total area in figure 7. The trends are clear: the lowest 
possible drag is obtained with low ratios of tail area to total 
area. 


CONCLUSIONS 


The influences of c.g. position and of tail plane design on 
the performance of a standard class sailplane have been shown. 

One important result is that the optimum c.g-. position is for a 
wide range of lift coefficients within the normal c.g. range. The 
calculations for the tail plane design have shown that the 
reduction of induced drag due to higher tail span is less 
important than the influence of profile drag. Low tail area and 
relatively low tail aspect ratio will give some advantages. It 
is remarkable that the best standard class sailplanes of today 
have a tail area of S^. ~ 1 m 2 and a tail aspect ratio of AR^. = 5 
while older sailplanes have for example S^_ = 1.5 m 2 and AR^_ = 6. 
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TABLE 1.- DATA SET OF A TYPICAL SAILPLANE 
[Reference loj 

Mass, kg 3 00 

Wing : 

Span, m 15 

2 

Area, m 10 

Mean aerodynamic chord, m 0.67 

Aspect ratio 22.5 

Pitching moment coefficient, zero lift -0.1 

Horizontal tail: 

Span, m 2.4 

Area, m 2 1.0 

Aspect ratio 5.76 

Distance between wing and tail aerodynamic centers, m 3.85 

Stability margin -0.15 
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DRAG COEFFICIENT 
















TAIL SPAN b t [m] 



a 1.5 2 3 4 

TAIL SPAN b*[m] 


Fig. 5: Variation of tail span, tail area = const. 
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